p-Wave Superconductivity near a transverse saturation field 
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We investigate reentrant superconductivity in an Ising ferromagnetic superconductor URhGe 
under a transverse magnetic field h x . The superconducting transition temperature for p-wave order 
parameters T sc is calculated and shows two domes as a function of h x . We find strong enhancement 
of T sc in the high-field dome near a saturation field h s where the spins align in the transverse 
direction. Soft magnons generate strong attractive interactions there. Spin components of the 
pairing show a significant change between h x <h B and h x >h s . We also discuss the appearance of 
superconductivity with zero-spin pair due to cancellation between external and exchange fields. 

PACS numbers: 74.20.-z, 74.20.Mn, 74.25.-q 



Ferromagnetic superconductivity (SC) has attracted 
much attention in condensed matter physics in the 
last decade^ U-based heavy-fermion compounds such as 
UGe 2 , URhGe, Ulr, and UCoGe show unconventional SC 
within their ferromagnetic phases^— Non-unitary SC's£ 
are believed to appear in these compounds and both the 
SC and ferromagnetism are caused in the / electrons at 
U sites. Their pairing mechanism and symmetry as well 
as novel self-induced vortex states are central issues to 
be clarified in the modern theory of unconventional su- 
perconductors. 

Two isomorphic compounds, UTGe(T=Rh, Co), ex- 
hibit SC at ambient pressure within their ferromag- 
netic state and have a similar Ising type anisotropy 
of magnetization*^ Spontaneous moment appears par- 
allel to the c-axis and its magnetization curve ex- 
hibits meta-magnetism when magnetic field H is ap- 
plied to one of the hard axes (b axis) with a notable 
mass enhancement £ This meta-magnetism is particu- 
larly prominent in URhGe and the moment gradually 
tilts with field and finally aligns parallel to H at 12 T£ 
Superconductivity appears below the transition tempera- 
ture T sc =0.24 K for H=0, and T sc decreases with H and 
disappears at 2 T for H || b£ Interestingly, SC reappears 
above 8 T and shows the highest T sc =0.42 K at 12 T£ 

Theories of ferromagnetic SC have been developed by 
various groups^ - — especially for UGe2*^ r — In this pa- 
per, we focus on the reentrant SC in URhGe, which has 
not been investigated theoretically. The physics of the 
reentrant SC is quite different from that in UGe2^~— 
The point is the presence of soft magnons in the Ising 
systems with transverse fields. We clarify it by extending 
the Scharnberg-Klemm (SK) theorj*2£ in the presence of 
ferromagnetism and transverse magnetic fields. In partic- 
ular, we demonstrate (i) strong enhancement of T sc near 
the saturation field, and (ii) strong field dependence in 
spin components (d- vector) of p-wave pairing. 

We consider a system of conduction electrons as shown 
in Fig. [TJ coupled to ferromagnetic localized spins, which 
have a magnetization M in the z-direction (correspond- 
ing to the c axis in URhGe at H=0). We apply a trans- 
verse magnetic field H x along the b axis and study its ef- 
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FIG. 1: (Color online) Schematic picture of the model. 

fects on ferromagnetic SC. For conduction electrons, we 
use isotropic dispersion £ p =p 2 /2m*, ignoring complex 
band structures^ Here, p and m* are the momentum 
and the effective mass, respectively. 

To model spin fluctuations in URhGe, we employ a 
simple model that can describe their essential behavior, 
and that is the ferromagnetic XXZ model.— To discuss 
SC, we will later introduce couplings between localized 
spins and conduction electrons. We treat local and itiner- 
ant degrees of freedom separately and this is qualitatively 
justified in the "duality" pictured There, ordered mo- 
ments are due to the incoherent part of / electrons, while 
itinerant physics such as SC is described by quasiparticlcs 
with heavy effective mass m*. In this paper, we denote 
magnetic moments and indices of the Kramers pairs sim- 
ply as "spin," but remember that these are complex com- 
binations of orbital and spin moments in heavy-fermion 
systems. 

The ferromagnetic XXZ model for localized moments 
is given as 
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where the spins Si with Sf=S(S + 1) are defined 
on a three-dimensional cubic lattice and h x =\xH x j S . 
fj,=0Afj,B is the value of magnetic moment measured in 
experiment^ where /is is the Bohr magneton. We have 
introduced spatial anisotropy Jfj as well as Ising spin 
anisotropy Jfj>J±>0 as shown in Fig. [TJ since the crys- 
tal structure can be regarded as coupled-chains running 
in the y(a) direction: J*,=J Z for bonds along x or z, 
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while rJ z along y with r > 0. This spatial anisotropy 
lifts degeneracy in p-wave superconducting order param- 
eters as will be discussed later. Although the tetragonal 
symmetry of the exchange couplings is higher than the 
orthorhombic one in URhGe, this is not essential for our 
discussions. 

We analyze the model ([1]) for zero temperature T=0 
by the linearized spin-wave approximation^ Due to the 
Ising anisotropy, the spontaneous moment points along z 
at h x —0. As h x increases, the magnetization tilts toward 
x with polar angle 6=sm~ 1 (h x /h s ) up to the saturation 
field h s =(4 + 2r) J z — 6J±. Magnons defined below de- 
scribe transverse fluctuations around this tilted direction. 

The magnon energy dispersion is expressed as 

4i)2 where eq=(cj-sl/2)h s +h x s g - J±(cg 
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with 



+I)7q-^*e7& 2« q =«g(/»./2+J,7«[-Jj.7 q ) 
7q=E/= a; , y , 2 ( C0S ffi- 1 ) 5 7^=7q+(^-l)(cosq. y -l), 
and ce(se)=cos9(s'm8). The dispersion has an en- 
ergy gap Syjh 2 - h 2 x for h x <h s and S^h x (h x - h s ) 
for h x >h s . At h x =h s , the excitations are gapless 



projection of Cooper pairs along h c . Different LL's for 
the p x state are decoupled to each other, while they are 
coupled for the (p y ,p z ) components. In the SK theory, 
the linearized gap equations are given as2£ 
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Eqr^o^S^J h s J z (q 2 + rq'y + q 2 ). Note that the model 

([T]) is ferromagnetic, and thus the linearized spin-wave 
approximation is fairly good even near h s . Since the 
magnitude of S does not play an important role in our 
approximation, wc set S=l. The renormalization of 
moment due to the quantum fluctuation is small and we 
neglect it hereafter. 

Let us investigate the interactions of magnons and 
conduction electrons (quasiparticles) . The quasipar- 
ticles feel the external field and exchange field via 
anisotropic antiferromagnetic Kondo coupling K XtVtZ : 
J2iv-Kv s iSi with sj' being the spin of the quasipar- 
ticle at the site i. Again, note that this "spin" 
is, indeed, "pseudospin" . Thus, in the mean field 
level, their energy dispersion is £ p ±=e p ± |h c |/2 with 
h c =(h*, 0, h*)=(n c H x -K x s g ,0, -K z c e ). Here, \i c =g* [i B 
with the effective g-factor of quasiparticles g* . 

The lowest-order fluctuations beyond the mean- 
field approximation yield a magnon-quasiparticlc 
interaction: X)ii/ w f a I*i'+-H-C-, where a\ is a magnon 
creation operator at the site i and s\ is the quasi- 
particle spin rotated such that the new z-axis is 
parallel to — h c . The coupling constant v v is given as 
(v x ,Vy,v z ) = (K z seh*+K x ceh z c ,-iKy,K z s g h z c -K x c g h*)/^ 
with h x c ' z =hl' z /\\x c \. 

Next, we derive linearized gap equations for p-wave 
SC— - and determine the transition temperature T sc and 
upper critical field H c2 - Since we are mainly interested 
in high-field states, we neglect small 47rM less than ~ 
1 kG in the magnetic induction^ This must be taken 
into account near h x =0 and is important for discussions 
about self-induced vortex states, which is one of our fu- 
ture problems. Since H || x, the gap equation for the 
p x state A a jv=[A^y, A xff , A xN ] T and the other orbitals 
{A y jv, A z tv} are decoupled, where N represents the Lan- 
dau level (LL), and "++" , " " , and "0" denote the spin 
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Here, V a is an interaction matrix of the spin part [see 
Eq-©] for the p a state. S includes both the Zeeman 
and orbital pair breaking effects and is diagonal in spin 
space. 

V a is determined from the one-magnon exchange pro- 
cesses and we calculate this with the static approxima- 
tion for the magnon Green's function D q (iw m )=l/(zu; m — 
Eq)^— 1/-Eq. Averaging all the momentum dependence 

assuming p-wave gap functions 77p" , V a is given as 
1 



V a =A±v 2 z \ w 2 




= A+v 2 z V+H+A-v 2 v V., 
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(3) 



where w=v x /v z and the first and the second columns 
of V a represent the equal spin pairing (ESP) with the 

spin function \-\ — h) and | ), respectively, while the 

third one represents the zero-spin pairing (ZSP) state 
(|H — ) + | — h))/\/2. Note that "-" and "+" denote spins 
parallel and antiparallel to h c , respectively. Here, 
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where the pairing form factor is chosen as rip=y/2 sinp a 

and this corresponds to V3p a in continuum systems.— 
p _k is related to the Bogoliubov transformation: 
tanh 2</>p_k=2fp_k/ep-k, and Nl is the number of sites. 

Before showing numerical results, we analyze V a for 
h x ~ h s . Since i? q is small around q=0, T sc is expected 
to be enhanced. It is important to note that the direction 
of the local moments continuously changes toward tt/2 for 
h x < h s , while 8~tt/2 above h s . Expanded up to first 
order in 56~ir/2 — 9, V a is given as 







A^K 2 ~ A a K 2 
A+K z — A a K y 

-rA+66 TAtSO A+K 2 + A-K% 



-TA+S9 

TA+80 | ,(5) 



where T=y/2K 2 z {[^ c h s / \ [ iK z ) - K./K,}- 1 - K x /K z }. 
For h x >h s , 66=0, and the gap equations for ESP and 
ZSP are decoupled. Note that the interaction is strongest 
for ZSP, since A^>0- This is natural, since the spin fluc- 
tuations are transverse, which scatter quasiparticles with 
± spin into =F spin states. 

Low-field longitudinal fluctuations^ also mediate in- 
teractions and we introduce a phenomcnological fer- 
romagnetic Ising interaction between quasiparticles: 
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- Jcc s i s j' witri J?c>® an d «e{x, y, z} is the direc- 

tion of ij bond. As before, we include spatial anisotropy: 
J l cc=Jcc for bond || x or z, while r'J cc for bond || y. 
Then, the interaction kernel © is replaced as V a —> 
V a + (h c z f J^ c V + (h c x /h c z )/4:, with keeping the form © es- 
sentially unchanged. 

Let us discuss the phase diagram of SC in the T-h x 
space determined by our calculations. Figures [2] and [3] 
show the /^-dependence of the SC transition tempera- 
ture of the p x - and p yjZ -states for the isotropic (r=r'=l) 
and anisotropic (r=1.2,r'=2) cases, respectively. Figure 
0]shows the spin part of the pairing, and will be discussed 
later. Our calculations show two domes of T sc in the 
phase diagram, and this qualitatively reproduces the ex- 
perimental datai2 One important observation is a strong 
enhancement and the presence of singularity in the high 
field dome at h x =h s . This manifests the change in the 
pairing state. The variations in the spin components are 
due to competition between the interaction strength and 
the Pauli depairing effects (PDE's) and we will examine 
this point in detail in the following. 

Let us investigate how the PDE's influence pairing 
states. The PDE in the SK theory originates in the ma- 
trix S in Eq.@. For h x >h s , the ESP's and ZSP are 
decoupled as discussed above, which simplifies analysis. 
The PDE is weak when the exchange field nearly cancels 
the external field h x c ~ 0. This is realized at high field, 
and we expect there the ZSP state, since the correspond- 
ing interactions are strongest. This is the Jaccarino-Peter 
effect,— which was originally proposed for a rare earth 
ferromagnetic metal. When the cancellation is not suffi- 
cient and thus PDE is strong, the ZSP state is suppressed 
and the ESP state is favored. If the interactions are much 
stronger than the effective magnetic field, two ESP com- 
ponents |H — h) and | ), have nearly equal amplitude. If 

the interactions arc not so strong, the ESP state with 
larger density of states (DOS) dominates 

For h x < h s , the situation is more complex. When 
PDE is strong, the ESP with large DOS is realized as 
expected also for low field, while for weak PDE all the 
spin components contribute to the superconducting con- 
densation, since the offdiagonal elements in V a are finite. 

Now, we explain the details of Figs. [2HH Our choice 
of parameters are: a mean-field ferromagnetic transi- 
tion temperature T r MF = J z (2 + r) = f K, H s =h s /n=W T, 
m*/m=40, r'J cc =283 K, K Z =W9 K, and K XtV =0.3K z , 
where m is the bare electron mass. Electron density is set 
to 1 per unit cell and the lattice constant^ is set to 3.5 A. 
Using these parameters, T sc at h x =0 is ~0.6 K for both 
parameter sets in Figs. [5] and g* is the control pa- 
rameter in our analysis below, which varies the Zeeman 
energy for quasiparticles. 

Let us first examine the orbital part of pairing. The 
orbital motion couples only to the external field H x as 
we mentioned above. Therefore, the p x -pairing is always 
decoupled from the other orbitals and different LL's are 
decoupled there, while the p y and p z orbitals are mu- 
tually coupled with LL's also coupled. In the isotropic 







J-sc 




T yz 




r =r'= 1 h 




< g* 2 1 





1 



FIG. 2: (Color online) Dependence of transition temperature 
Tgc yz on transverse field h x for r=r'=l and g*—2. 




h x /h s 

FIG. 3: (Color online) ^-dependence of for the 

anisotropic case r—1.2, r'=2 and five sets of g* . The change 
of low-field dome with <?*'s is negligible. For h x >h s , the ZSP 
state is decoupled from the ESP state. Its Tf* is shown by 
lines without symbols, except for g* = l,2, and 6, where T]f* 
is nearly zero. 

systems, the polar p^-state with A^=0 has the highest 
H C 2 as shown in Fig. [2] We note that the N=Q LL plays 
a central role in both p x and p y . z states^ 

When the interactions arc spatially anisotropic in 
space with r, r'>l, the p x state is suppressed and a p y - 
dominant state becomes stable, and vice versa for r, r'<l. 
In the following, we analyze the former case, which is rele- 
vant for URhGe. Figure [3] shows T»* for r=1.2 and r'=2, 
as a typical case of strong anisotropics, for several values 
of g*. Here, Tf c <C Tj^ for all the parameters and not 
shown. For the low-field SC, a p y dominant state has a 
lower energy even for small r' > 1.05, while less sensitive 
to r. This is because that the magnon gap is sufficiently 
large there and thus the transverse fluctuations are neg- 
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FIG. 4: (Color online) Spin part of p y pairing in the iV=0 
LL for r=1.2 and r'=2: A++(0), A~-(A), and A£ ( V ) 
at T=TJ / C f. Lines represent |h c |. (a) Low-field part, g*- 
dependence is negligible, (b)-(f) High-field part near h s for 
various g*'s. Localized moment (long arrow) and quasiparti- 
cle spin (short one with dot) are schematically shown. The 
magnitude of the quasiparticle spin is small < 0.05 near h s . 

ligiblc. For the high-field SC, both r and r' affect its T sc . 
Increasing r enhances T^*, while it suppresses Tf c . For 
the dependence on r' , increasing r' suppresses both T sc 's. 
The suppression is stronger for Tf c than for T)j* because 
the attractive force does not decrease for p y pairs. It is 
also possible to realize a transition in the SC from p x 
to p V:Z symmetry or vice versa, as h x increases, which is 
determined by the details of the anisotropies. 

Second, we discuss the spin part of pairing, which is 
related to the d- vector of SC, and examine for the data 
in Fig. [3] the Zeeman effects with controlling g*. Figure 
2] shows the change in the spin part of pairing with h x for 
the orbital component p y and N=0 LL. They are calcu- 
lated just at TJf* and normalized as z^atScJ Aaiv| 2 =l- 
Figure H^a) shows the change in the low-field dome. We 
find that the g*-depcndence is negligible and the pair- 
ing is the ESP state with only one type of spin, which is 

parallel to the internal field h c , | ). This is consistent 

with the proposed gap symmetry for URhGe^ Figures 
IUb)-[lIf ) show the results in the high-field dome. The re- 
gion of h x >h s is, in most of the cases, also the pure ESP 
state, but now the two spin components have almost the 
same amplitude. This is because |++) pairs are scat- 



tered only to | ) pairs and vice versa, when h x >h s . 

However, their amplitude difference may become larger 
depending on g* and |h c |. In the region of h x <h Sl a 
small amplitude of the ZSP component hybridizes. The 
two ESP components are still dominant, but their differ- 
ence is now larger and the /i^-dependencc is noticeable, 
particularly for larger g*'s. For all g*'s, transverse spin 
fluctuations near h s favor the linear combinations of two 
ESP states, while ZSP is suppressed by the PDE, except 
for some regions in Figs. Hfd) and HJe) as we discuss 
below. 

The quasiparticle spins respond to both the exchange 
field of the localized moments and the external magnetic 
field, and g* controls the latter part. Thus, varying g* 
affects the angle between the quasiparticle and the lo- 
calized spins as schematically depicted in Figs. lUJb)- 
EJf). When the external and exchange fields are canceled, 
h c =0, the quasiparticle spins feel no field. Indeed, the 
cancellation can become nearly perfect for g*=4.2 and 5. 
This leads to emergence of ZSP states (Jaccarino-Peter 
effect) and occurs in the shaded regions in Figs. 0Jd) 
andQIc). When h s is close to the point |h c |~0 as in Fig. 
|4^e), the expansion in Eq.Q is not sufficient, and this 
results in the almost vertical slope of A y o for h x < h s , 
but this situation is not of main interest in this paper. 

Let us finally comment about URhGe. We have suc- 
ceeded in reproducing the overall phase diagram and two 
SC domes as shown in Figs. [2] and [31 The transition tem- 
perature for the second dome shows strong enhancement 
near h s . The superconducting gap symmetry is basically 
equal-spin p^^-pairing for both domes. It is | ) for 

the low-field, and - (|++) + | ))/v / 2 for the high-field 

SC. The gap functions have a point node due to the finite 
ip z component, but since its amplitude is very small, it 
might be difficult to experimentally distinguish this from 
that with a line node. Detecting the change in the spin 
components of pairing near h s is an experimental test for 
the present theory. Exploring the sudden change in pair- 
ing symmetry due to the Jaccarino-Peter effect is another 
interesting challenge in this field. 

For a more precise theoretical analysis, one should take 
into account the optimization of form factor ?y£, longitu- 
dinal spin fluctuations for h x < h s ^ the weak first-order 
transition^ the variations of the Fermi surfaces near h s ^ 
and anisotropy in the effective mass^ which might im- 
prove quantitative agreement between our theory and ex- 
periments. Our main conclusion for the superconducting 
order parameters for h x ~h s is that the pairing interac- 
tion prefers mixing |H — h) and | ) and this is in clear 

contrast to nearly pure | ) state, which is selected by 

Zeeman energy. Difference in their amplitudes depends 
on the details, and a more quantitative analysis about 
this is one of our future plans. 
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